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Abstract 

The two-loop contribution to the electric dipole moment (EDM) and the chromo 
electric dipole moment (CEDM) of an arbitrary fermion / induced by the most gen- 
eral renormalizable tbW coupling with complex left- and right-handed components 
(ol and a^) is calculated. The analytical expressions are numerically evaluated and 
the current experimental constraints on the electron, neutron and mercury atom 
EDMs are used to obtain a bound on the complex phase Im(a2a_R)- It is found that 
the most stringent constraint, lui{a\aii) < 2.33 x 10~^, arises from the neutron 
EDM. 



Although there is experimental evidence of CP violation, its origin still re- 
mains a mystery. In the standard model (SM), the only source of CP violation 
is the Cabbibo-Kobayashi-Maskawa (CKM) phase, which appears to be the 
origin of the CP violating phenomena observed in nondiagonal processes in- 
volving the K and B mesons On the other hand, as diverse studies had 
suggested [2], the CKM phase has a rather marginal impact on flavor- diagonal 
processes, such as the electric dipole moment (EDM) of elementary particles. 
For instance, the EDM of fermions arises up to three loops within the SM 
thereby being extremely suppressed. It can be significantly enhanced, how- 
ever, in several SM theoretical extensions, in which it can be induced at lower 
orders via new sources of CP violation. It means that any experimental signal 
associated with an EDM would point to new physics. Therefore, the EDM of 
light fermions, such as the electron and neutron, has been the subject of con- 
siderable interest in theories beyond the SM, such as supersymmetric models 
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[1], multi-Higgs models [5], left-right symmetric models [6], and other theories 
[7j. Along these lines, a potential source of CP violation may be the tbW cou- 
pling, whose study will be a top priority at the CERN large hadron collider 
(LHC). Such CP- violating effects would be induced via a complex phase aris- 
ing from the simultaneous presence of both left- and right-handed components 
in the tbW vertex, a scenario which is predicted indeed in several SM exten- 
sions. The purpose of this work is to calculate the two-loop contribution of 
the tbW coupling to the EDM of a fermion. Although we will obtain a result 
valid for any charged fermion, our main goal is to use the experimental limits 
on the EDM of the electron and the neutron to constrain the complex phase 
associated with the anomalous tbW vertex. As a byproduct, we will obtain the 
chromo electric dipole moment (CEDM) of / and use the current experimental 
bound on the EDM of the mercury atom to constrain the anomalous part of 
the tbW vertex. 

The most general renormalizable tbW coupling is given by the following La- 
grangian 

C = -j=t{aLPL + aRPn)bW+ + H.c, (1) 



where Pl,r are the usual left- and right-handed projectors and are un- 
known complex coefficients. The contribution to the on-shell ff'j vertex, which 
defines the electromagnetic properties of the fermion /, is given, in the uni- 
tary gauge, via the Feynman diagrams shown in Fig. [H Similar diagrams, 
with the photon replaced by a gluon, give rise to the CEDM. We will present 
the calculation of the ffj vertex and then generalize the result to the ffg 
coupling. 

From the amplitude for the on-shell ffj vertex we can obtain the EDM of 
/, which is the term proportional to the Lorentz tensor structure i'j^cr ^uQ'^ ■ 
We will use the Feynman parameter technique to integrate over the arbitrary 
internal momenta. For the sake of completeness, we will present the most 
relevant details of the two-loop calculation in appendix A. 

After solving the four- dimensional integrals arising from Fig. [H the EDM of 
/ can be expressed as follows 

df = - (^] (£-) Ar,x,x,xjm(ala«)/, (2) 



where Xa = ma/mw, A''^ = 3 is a quark color factor, and / stands for the 
quadruple integral given in appendix A [Eq. (1A.25P ]. which must be numeri- 
cally evaluated in the most general scenario, namely, for internal and external 
fermions with non- negligible masses. On the other hand, in the Xf = Xi = 



2 



/(Pl) /. /(P2) /(Pi) !i f(pi) 

Fig. 1. Two-loop diagrams contributing to the on-shell //7 vertex. There are other 
diagrams in which the photon is emitted from the W boson or the internal fermion 
/i, but they do not contribute to the EDM of /. The contribution to the CEDM 
arises from a similar diagram with the photon replaced by a gluon. 

approximation, which is suited for the purpose of our work as we are inter- 
ested in evaluating the EDM of very light fermions, the parametric integral 
I adopts the more simple form given in Eq. (1A.27I) . which can be simplified 
further after some algebra: 

J Qb + Qt , Qtxl-Qbxl ^ (^t\,}^ wi\ 



with 



Qb x{l-x) 
2 (gh - x{l - x)y 



log 



9b 



x{l 



x] 



(4) 



It is not possible to integrate the F^lx) function in terms of elementary func- 
tions, but it is straightforward to obtain a numerical solution due to the fact 
that F^lx) is well-behaved in the corresponding domain. 

It is useful to express the numerical value of the EDM and the the CEDM for 
a light fermion. Inserting the approximate value of the integral /, we obtain 
for a very light fermion: 

\df \ = 3.194 X 10~^^ Xflmialan) e ■ cm, (5) 



As for the CEDM of /, dj, it follows easily from the above results. We just 
need to make the following replacements df = df{e — > gs,Qt ^,Qb ^ !)• 
Numerical evaluation gives: 

\df \ = 3.056 X 10~^^x/Im(a2ai?) gs ■ cm, (6) 



We turn to analyze our results in the light of the current experimental limits 
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on the EDM of fermions. Let df{Exp) be the experimental hmit of the EDM of 
a fermion /, which is not necessarily a light one. Then, our theoretical result 
for df can be translated into the following bound on the complex phase: 

a„) < ( *-^] C-!^] ( ( l''^<fP>h . (7) 

\ o;^ / V e / \mbmtmf J \ I J 



This expression is exact as long as Eq. flA.27p is used for the parametric integral 



/. In the case of a light fermion, we can use the result of Eq. ([5]). 

We are now ready to constrain the phase lm{a*j^af>). For this purpose, we will 
use the experimental data on the EDM of the electron, the neutron and the 
mercury atom. 

The current experimental limit on the EDM of the electron is [Sj : 

|4(Exp)| < 7 X 10-2® e • cm, (8) 

Using this limit, we obtain the following bound 

lm{alaR) < 0.345. (9) 

As far as the neutron is concerned, there are three different approaches to 
estimate its EDM dn [2]- In the chiral Lagrangian approach, dn is expressed 
in terms of the quark CEDMs, whereas in the QCD sum rule approach it is 
expressed as a combination of both quark EDMs and CEDMs. Due to the 
large discrepancies arising from these two approaches, it is convenient for the 
purpose of this work to estimate the neutron EDM using the non-relativistic 
SU (6) quark model along with naive dimensional analysis for the quark CEDM 
contributions. In this approach, we have 

dn = ^{4:dd-du), (10) 



where 



6 ~ 

du,d = Vdu^ + V^du^d (11) 



with 0.61) and r7(~ 3.4) being the respective QCD correction factor from 
renormalization group evolution. The CEDM contributions should be included 
as long as they are of similar size to those from the EDMs, otherwise they can 
be neglected. In our case, the quark CEDM is of similar size of the EDM, but 
its contribution to Eq. (ITTi) is suppressed by more of one order of magnitude 
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due to the factor l/(47r). It is thus safe to neglect the CEDM contribution. 
On the other hand, currently the most stringent bound on the neutron EDM 
is [9j: 

\dn{Exp)\ < 2.9 X 10-2^ e ■ cm. (12) 



As usual, we take rriu ~ rrid ~ rnn/<i, with m„ the neutron mass. The above 
Eqs. leads to a stronger bound than the one found from the electron EDM: 

lm{alaR) < 2.33 x 10^1 (13) 



As for constrains on the EDM of diamagnetic atoms, the most stringent one 
can be obtained from the mercury atom. A constraint on the Im(a2ai?) can be 
found following the approach of Ref . [lOj , in which the experimental constraint 
on the mercury EDM is translated into the bound \du — dd\ < 2 x 10~^^ 
cm. Explicit calculation shows however that the resulting upper constraint on 
Im(a2aij) is weaker than the one obtained from the electron EDM by about 
one order of magnitude. So we will not consider this constraint in this work. 
We can conclude that the most stringent bound is the one obtained from the 
neutron experimental data. 

We now would like to compare our constraint with other ones appearing in 
the literature. It has been customary to express the left- and right-handed 
parameters in the following way: 



= 1 + /^Le^'^^ (14) 
aR = KRe'^^ (15) 

where k^ ji and 4'l,r are real parameters. Moreover, k^ ji > 0. This parametriza- 
tion of the tbW vertex is simply the SM contribution [ql = 1, a_R = 0) plus an 
anomalous complex term expressed in polar form. It follows that our constraint 
(fT3ll translates into 

kr sin (pR + klKr sin(0H - (pi) < 2.33 x 10"^. (16) 



Constraints on these parameters have already been reported in the literature. 
Data from B meson physics allowed the authors of Ref. [TT] to impose the 
following limits: 



kl sin (f)i<3 X 10~ 



(17) 
(18) 
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On the other hand, the CLEO Collaboration data on the decay 6 — 57 have 
been used in Ref. [12] to constrain the right-handed parameters: 



0.0035 < KRCOs(f)R + 20k% < 0.0039 x 10^^ (19) 



However, as pointed out in Ref. [12], this constraint is not sensitive to CP- 
violating effects, which are the ones we are interested in. 

It is worth combining the above constraints to find the allowed region on 
the (pL — (pR plane. Since Eq. (fT3l) depends on four unknown parameters, we 
will assume appropriate values for hl,r- In particular, the bounds hl < 0.01 
and kr < 0.2 have been derived from the b ^ sj decay [I3]. These bounds 
are somewhat restrictive and thus larger values of these parameters are still 
possible. It is thus interesting to analyze the following scenarios: {i)nL ~ i^r 
and (ii)KL ^ kr. In the ^ i^r scenario our constraint is redundant as (1181) 
gives a tighter constraint on kr sin (pR. In Fig. 2 and 3 we have plotted the area 
on the (pL — (pR plane allowed by the constraints f|T3l) and f|T7|) - f|T9l) for various 
values of and hr. The surviving region obtained after combining all the 
constraints is also shown. It is easy to see that the CLEO constraint is very 
constraining on the qr parameter, although it is not useful to constrain the 
kl parameter, let alone the (pi phase. The main advantage of our constraint is 
that it is useful to constraint the CP violating phase and it is expected to give 
a very stringent constraint on this phase once the experimental constraint on 
the neutron EDM is improved. 




Fig. 2. Allowed area on the (pL — (pR plane as obtained from the constraints (jl3p and 
(|17|) -(|19|). The unshaded area represents the area allowed by the constraint from 
the neutron EDM (upper left), B meson physics (upper right), the CLEO data on 
6 — > 57 (lower left), as well as the surviving region obtained by combining all the 
constraints (lower right). The values = 1 and kr = 0.05 have been used. 
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Fig. 3. The same as in Fig. 2 for = kr = 0.05. 

We would like to emphasize some advantages of our method for constraining 
the tbW coupling. First of all, in obtaining our constraint no extra assumptions 
were made. In fact, although the result arises from a two- loop calculation, it 
is free of ultraviolet divergences and the expression given in Eq. ([7j) is exact 
since the parametric integral / can be numerically evaluated for any external 
fermion, including even a very heavy one. Our result thus can be useful for 
predicting the EDM of the /i and r leptons or the s and c quarks, for instance. 
In addition, our bound can be easily updated. Along these lines, in Ref. [H] 
a proposal was presented for improving the experimental limit of the EDM 
of the electron by 3 orders of magnitude. This would lead to a bound on 
Im(a2a_R) of the order of 10"'^, which is of the same order of magnitude than 
that arising from B meson physics. More recently, the nEDM Collaboration 
[15j has presented a proposal [16j to improve the current limit on the neutron 
EDM by 2 orders of magnitude. This potential constraint would lower the 
bound on lm{a*^aii) to the level of 10~^, which would be more stringent by 
about one order of magnitude than those obtained from B meson physics and 
the CLEO data. 

In closing we would like to emphasize the relevance of the present work. Im- 
portant information on the origin of CP violation may be extracted from 
the measurement of the EDM of elementary particles. Several sources of CP 
violation are predicted in beyond-the-SM models, and the most stringent ex- 
perimental limits imposed on the EDM of the electron or the neutron would 
allow us to asses their relative importance or eventually to rule them out. In 
this paper, we have studied the impact of a complex phase associated with the 
most general renormalizable tbW coupling with both left- and right-handed 
components on the EDM of an arbitrary / fermion. For the sake of complete- 
ness, the calculation of the two- loop amplitude was analyzed to some extent. 
The resultant expression can be straightforwardly used to predict the sensi- 



7 



tivity of the EDM of a light or heavy fermion to a complex phase appearing in 
the tbW coupling. In particular, our theoretical result was numerically eval- 
uated in the scenario of a light fermion and the outcome was combined with 
the experimental limits on the EDM of light fermions to obtain an inequality 
that can be easily updated to constrain the tbW coupling. Using the most 
recent experimental constraints on the EDM of the electron and the neutron, 
it was found that the latter gives the most stringent bound on the tbW com- 
plex phase, which is one order of magnitude less stringent than those obtained 
from B meson physics. We would like to emphasize however that our bound 
could be improved by about two orders of magnitude if the neutron EDM is 
measured with more precision at a near future, as recently proposed |15j . 
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A The two loop calculation 



Using the SM Feynman rules and the one induced by the Lagrangian ([T]), the 
amplitude for the //7 coupling can be written as: 

M=uip2)T^uip^yiq,\), (A.l) 



where is the two-loop vertex function, which is given by 




(A.2) 



with QtiQb) the electric charge of the t{b) quark in units of the positron charge, 
rrii the mass of the internal fermion that couples to the W boson and the / 
fermion, and 



At = [kl - ml] + ^2 - Pif - m]] [{h + ^2 - - m% (A.3) 
A, = [kl - ml] [(fci + k2-vi?- ml] [(fci + k2- V2? - mH (A.4) 




,b 



Tr\^0{aLPL + ai?PR))(^2 + mfc)7„(a*Pi + a^P^^) x 
(Vi +^2 -l/i + mth^il^i +^2 -1/2 + mt)] , 
Tr[ff3(alPL + a*i^PR))(J^2 + mtha{aLPL + anPn) x 
(Vi +^2 -l/i + ^6)7^(^1 +^2 -^/2 + . 



(A.7) 



(A.8) 



This amplitude generates contributions to all the form factors associated with 
the on-shell //7 vertex, but we are only interested in the term proportional 
to ^5afj^u<f . The latter can be isolated after the identity ep^uaiil^l^ = —2750"^!/ 
is used. After some algebra, the contribution to the EDM of / can be written 
as 



where sw is the sine of the weak angle. Note that the integrals in D dimensions 
were written in four dimensions as they are free of ultraviolet divergences. 

We will proceed to solve the two-loop amplitude. Using Feynman parameters 
for the integral over k2, one obtains 




d% ( Qt Qb \ 



(A.IO) 



where 



Rt = x{x 
Rb — x{x 



r2 
t ) 



r2 
b ' 



(A.ll) 
(A.12) 



with x — xi + X2 and 




(A.13) 



X X 




(A.14) 
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M^^m^ + ^-'^. (A.15) 

Once the integral over k2 is done, Feynman parametrization for the integral 
over ki leads to 



rj = — —mbmt\m{alaR){'^5a^^q'')T{A) / dxi / 

1 l-yi l-y 



x{x — 1) 







where we have introduced the following notation: y — yi + y2, p — yip\ + 
y2P2 + /|/3, and 

M^, = M^^yg + m^y + m%y + y3)(y + - 1) + ^,^(1 -y- Vs)- (A.17) 



The EDM of / can thus be written as 

A /Tr-2 



where we have introduced the dimensionless variables Xa — nia/mw and 

1 l-xi 1 l-y / 

Qt _ Qb 



I ^ J dxi j dx2 j dyi j dy2 j ^^32^(1 - a:)(y + t/s) f ^ - ^ j , (A.19) 



0000 



with 



Ft = fty3 + f, (A.20) 
Fb = fm + I\ (A.21) 

ft = x{l-x)x} + {xl-xf)x + xl, (A.22) 

/5 = ,x(l - x)x} + (x^ - xl)x + x?, (A.23) 
/ = a;(l - x)y + x{l - x){y + y3){y + ys - l)xj 

+x{l - x){l - y - y3)xl (A.24) 

We can further simplify the parametric integral / . Explicit integration over 
ya leads to 

1 l-xi 1 1-t/i 

I = j dxi j dx2 j dyi j dy2Fi{x,y), (A.25) 
0000 
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where 



x(l — x) 

Jb 




-y){f - hy) 
f + fb{i-y) 




It is not possible to solve the integrals over yi nor y2 analytically, so numerical 
evaluation would be necessary. However, a relatively simple expression can be 
obtained in the xj = Xt = limit, which is a good approximation for the case of 
a light fermion doublet coupling to the W boson. This is suited for the purpose 
of our work as we will use the experimental limits on the EDM of the electron 
and the neutron to constrain the tbW coupling. In such an approximation, 
once the integral over yi and ?/2 are done, the parametric integral / reduces to 



l—Xl 



dxi / dx2F2(x) 



(A.27) 



where 



F2(x: 



^Qb 

2 



x{l — x) 
{gb - x(l - x)y 



log 



Qb 

x(l — x) 



gb - x{l - x) 



(A.28) 



with gh = fbixf = 0) = x^ + {xl — x^)x. The second term between the 
square brackets can be integrated straightforwardly, whereas the term with 
the logarithm cannot be integrated in terms of elementary functions. The 
final result is given in Eq. ([3]). 
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